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Denoting P{G, q) as the chromatic polynomial for coloring an n- vertex 
graph G with q colors, and considering the limiting function W{{G},q) = 
liuin^ao P{G,qy^" , a fundamental question in graph theory is the following: 
is Wr{{G}, q) = q~^W{{G}, q) analytic or not at the origin of the 1/q plane? 
(where the complex generalization of q is assumed). This question is also rel- 
evant in statistical mechanics because W{{G},q) = exp{So/kB), where 5o is 
the ground state entropy of the g-state Potts antiferromagnet on the lattice 
graph {G}, and the analyticity of Wr{{G},q) at 1/g = is necessary for the 
large-g series expansions of Wr{{G},q). Although Wr is analytic at 1/q = 
for many {G}, there are some {G} for which it is not; for these, Wr has no 
large-g series expansion. It is important to understand the reason for this 
nonanalyticity. Here we give a general condition that determines whether or 
not a particular Wr{{G},q) is analytic at 1/q = and explains the nonana- 
lyticity where it occurs. We also construct infinite families of graphs with Wr 
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functions that are non-analytic at 1/q = and investigate the properties of 

these functions. Our results are consistent with the conjecture that a sufficient 

condition for Wr{{G},q) to be analytic at 1/g = is that {G} is a regular 

lattice graph A. (This is known not to be a necessary condition). 
05.20.-y, 64.60.C, 75.10.H 



I. INTRODUCTION 



The chromatic polynomial P{G, q) of an n- vertex graph G and the asymptotic limiting 
function 

W{{G},q)=\m.^P{G,qf'- (1.1) 

play important roles in both graph theory |l|]- and statistical mechanics [0]- 0. Here 
P(G, q) is defined as the number of ways of coloring the graph G with q colors such that 
no two adjacent vertices have the same color, and {G} denotes the limit as n ^ oo of the 
familiy of n-vertex graphs of type G. The connection with statistical mechanics is via the 
elementary equality P(G, q) = Z{G, q,T = 0)paf, where Z{G, q,T = 0)paf is the partition 
function of the zero-temperature g-state Potts antiferromagnet (AF) on the graph G, 
and the consequent equality (in the n — oo limit) W{{G},q) = exp{So{{G},q)/kB), where 
Sq{{G}, q) denotes the ground state entropy of the g-state Potts AF on {G} (typically a reg- 
ular lattice, {G} = A with some specified boundary conditions). Given the fact that P{G, q) 
is a polynomial, there is a natural generalization, which we assume here, of the variable q 
from integer to complex values. Since an obvious upper bound on P{G, q) describing the 
coloring of an n-vertex graph with q colors is P{G^ q) < g", and hence VTdG}, q) < q, it is 
natural to define a reduced function 

Wr{{G},q) = q-'W{{G},q) (1.2) 

A fundamental question in graph theory concerns whether Wr{{G},q) is analytic or not at 
the origin, 1/g = 0, of the z = 1/q plane. This question is important in both graph theory 
and statistical mechanics because a standard method for studying this function or equivalent 
reduced W functions is to calculate a large-g Taylor series expansion about the point 1/g = 
||T0| - [1T2| . However, if WrdG}, q) is nonanalytic at l/q = 0, then one cannot carry out such 
a Taylor series expansion in the usual manner. Indeed, we recently discussed an example, 
namely, the bipyramid graph i?„, for which Wr{{B},q) is not analytic at l/q = |l^ (see 
also 0). 

Clearly it is important to understand better the differences between the families of graphs 
that yield Wr{{G},q) functions analytic at 1/g = and those which produce Wr{{G},q) 
functions that are nonanalytic at 1/g = 0. In the present paper we shall address this problem. 
We shall give a general condition that determines whether or not a particular Wr{{G},q) 
is analytic at 1/g = 0. This explains the source of the nonanalyticity in the cases where it 
occurs. We also construct infinite families of graphs with Wr functions that are non-analytic 
at 1/g = 0. These serve as a very useful theoretical laboratory, and we study the properties 
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of the resultant Wr functions in some detail. A salient point is that none of the {G} that 
we construct with Wr{{G},q) that are nonanalytic at 1/g = is a regular lattice graph 
{G} = A. Thus, anticipating the later discussion in this paper, our work is consistent with 
the conjecture that (in the n ^ oo limit) a sufficient condition that Wr{{G},q) be analytic 
at 1/g = is that {G} = A is a regular lattice graph. (We know from our previous work |]13[ 
that this is not a necessary condition.) We state it as a conjecture since we are not aware of 
any proof of the analyticity at 1/g = of Wr{A, g) for a regular lattice A in the literature. 



Indeed, in Ref. [0, it was acknowledged that there was no general theory for the existence 
of the limit (|Ll|) and, hence also, in our notation, the reduced function Wr, even in the case 
of regular lattices. 

Before proceeding, it is necessary to clarify the definition of W{{G},q) for values of g 



that are not positive integers. As we discussed in Ref. |T^, for certain ranges of real g, 
P{G, g) can be negative, and, of course, when g is complex, so is P{G, g) in general. In these 
cases it may not be obvious, a priori, which of the n roots 

P(G',g)i/" = {|P(G,g)|i/V™/"} , r = 0,l,...,n-l (1.3) 

to choose in eq. (|1.1| ). Consider the function W{{G},q) defined via eq. (|1.1|) starting 
with g on the positive real axis where P{G, g) > 0, and consider the maximal region in the 
complex g plane that can be reached by analytic continuation of this function. We denote 
this region as -Ri. Clearly, the phase choice in (|1.3| ) for g G -Ri is that given by r = 0, 
namely P(G',g)i/" = |P(G, g)!^/". However, as we showed via exactly solved cases in Ref. 
p!3| , there are many families of graphs {G} for which the areas of analyticity of W{{G},q) 



include other regions not analytically connected to Ri, and in these regions, there is not, in 
general, any canonical choice of phase in (|1.3|) . 

A second subtlety in the definition of ({G}, g) concerns the fact that at certain special 
points gs, the following limits do not commute [|l3l (for any choice of r in eq. (|1.3|) ): 



lim lim P(G,g)^/" ^ lim lim P(G,g)i/" (1.4) 

One can maintain the analyticity of W{{G},q) at these special points g^ of P{G,q) by 
choosing the order of limits in the right-hand side of eq. ( p..4| ): 

W{{G},qs)D,,. ^ lim lim P(G,g)^/'^ (1.5) 

As indicated, we shall denote this definition as Dqn, where the subscript indicates the order 
of the limits. Although this definition maintains the analyticity of iy({G}, g) at the special 
points g^, it produces a function VrdG}, g) whose values at the points g^ differ significantly 
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from the values which one would get for P{G^qsY^^ with finite-n graphs G. The definition 
based on the opposite order of limits, 

W{{G}, qs)D^^ ^ lim hm P(G, qf''' (1.6) 

gives the expected results like W{G, qs) = for = 0, 1, and, for G 3 A, g = 2, as well as 
W((tri)n,q = 3) = 1 (where (tri)n denotes a triangular lattice with n sites and boundary- 
conditions that do not introduce frustration for g = 3), but yields a function iy({G}, q) with 
discontinuities at the set of points {qs}- In our results below, in order to avoid having to 
write special formulas for the points g^, we shall adopt the definition Dgn but at appropriate 
places will take note of the noncommutativity of limits ( p. .40 . 



II. CONSTRUCTION OF FAMILIES WITH Wri{G},q) NONANALYTIC AT 1/g = 

A. General Algorithm and Calculation of Chromatic Polynomial 

In general, as discussed in Ref. [0, for a given family {G}, the corresponding PF({G}, g), 
at least as defined via the order of limits ( |1.5| ), is an analytic function in certain regions of the 
complex g plane. These regions are separated from each other by curves (or lines) comprising 
the union of boundaries B. W{{G}, q) is nonanalytic on these boundaries. Clearly, B is the 
same for W{{G},q) and Wr{{G},q). Because P{G,q) is a polynomial with real (actually 
integer) coefficients, it follows that B is invariant under complex conjugation, i.e., 

B{q)=B{q*) (2.1) 

A basic question is whether, for a given family {G}, some portion of the boundary B extends 
to complex infinity in the g plane, so that Wr{{G},q) is nonanalytic at 1/g = in the 1/g 
plane. Related to this, an important question is whether there is a general algorithm for 
producing a family {G} of graphs such that in the n —>■ oo limit, the boundary B extends to 
complex infinity in the g plane. We answer this question in the affirmative and present the 
following algorithm. 

Consider a family of graphs {G}. If this family already has the property that the limiting 
function W^({G}, g) has a region boundary B that extends to complex infinity in the g plane 
(i.e., to 1/g = in the 1/g plane), then we have no work to do to get such a boundary. So 
assume that {G} is such that W{{G},q) has a region boundary B that does not extend to 
complex infinity in the g plane. As discussed in (section III and Theorem 1, eq. (3.1) of) 
Ref. [|13|, a rather general form for the chromatic polynomial of a graph G is 
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P(G„,g) = g(g- l){co(g) +^c,(g)a,(g)"} (2.2) 

i=i 

where Cj(g) and aj{q) are polynomials in q. Here the aj(g) and Cj-^o{q) are independent of 
n, while co{q) may contain n-dependent terms, such as (—1)", but does not grow with n 
like a". Obviously, the reality of P{G,q) for real q implies that Cj{q) and aj{q) are real for 
real q. The condition that B does not extend to infinite distance from the origin in the q 
plane is equivalent to the condition that for sufficiently large |g|, there is one leading term 
aj{q) in eq. (|2.2| ). Here we recall that "leading term a^(g) at a point q" was defined in 
Ref. |jl3| as a term satisfying |a^((3')| > 1 and |ci£(q')| > for j ^ i. (If the cq term is 

absent and Na = 1, then the sole ai(g) may be considered to be leading even if |aj(g)| < 1.) 
We require sufficiently large |g| so that, for our analysis, there is a switching between only 
two leading terms ae. In principle there might be such switching between more than two, 
so that B would include more than two components running to complex infinity in the q 
plane. However, for the families that we have constructed via our algorithm and studied, we 
find, for sufficiently large \Im{q)\, only two such components. As required by the symmetry 
( p.l|) , these components are mapped to each other under complex conjugation. Now adjoin a 
complete graph Kp to Gn in such a way that each vertex in Kp is adjacent, i.e., connected by 
bonds (=edges), to each of the vertices of On- Here, recall that a vertex graph is defined 
as "complete" and labelled Kp if each vertex is completely connected by bonds with all the 
other vertices of the Kp graph. Denote the resultant graph as {Kp x A basic theorem 
of graph theory states that if a graph H is obtained by adjoining a vertex to a graph G such 
that this point is adjacent to all of the vertices of G, then the chromatic polynomials are 
related according to 

P{H,q) = qP{G,q-l) (2.3) 

Applying this iteratively p times, we obtain the result that 

p-i 

P{Kp X G„, q)=[Y[{q-s)] P{Gn, q - p) (2.4) 

s=0 

Next, we select one vertex in Kp and remove b bonds connecting it to other vertices of 
Kp. Since each vertex of Kp has p — 1 bonds connecting it to other vertices of Kp, this implies 
that we can only remove this many such bonds, i.e., 

l<b<p-l (2.5) 

We denote the resultant graph as {Kp x Gn)rb, where the subscript signifies the above removal 
(r) of b bonds. In order for this to be nontrivial, i.e., for 6 > 1, we thus require that 
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p>2 

The conditions ( p.6| ) and ( p. 51) will be assumed henceforth. Using eq. 
of the addition-contraction theorem [|T5|, we obtain the important result 



(2.6) 

and r applications 



P{{K, X Gn)rt, q) = P{K, X G^, g) + bP{K,_, x G^, g) 

p-2 



[mi 

s=0 



- S] 



{q-{p- 1))P(G„, q-p)+ bP{Gn, q-{p- 1)) (2.7) 



This is our general formula for the chromatic polynomial of {Kp x Gn)rb, ioT an arbitrary 
n- vertex graph Gn- Substituting the expression (P^), we obtain 



P-2 



P((irpxG0.fe,g) = [n(9-^) 

r Na 

{q - p + l){q - p){q - p - l){co(g - + ^ ^(g - p)aj(g - p)"} 



Na 



+b{q -p+l){q- p){co{q - P + 1) + ^ c^(g - P + l)aj{q -p + 1)"} 



(2i 



B. Boundary B for {(i^p x G)rb} 
We denote the n — > oo limit of the families Kp x G„ and (Kp x Gn)rb as 



and 



lim Kp X Gn = {Kp x G} 



lim {Kp X G„)^,fe = {(iTp X G),.^} 



(2.9) 



(2.10) 



respectively. As discussed in Ref. [0, the boundary B for iy({G}, q) is the locus of points 
in the q plane where there is a switching between different leading terms a£ in eq. ( |2.2|) . 
Since {G} was assumed not to have B extending to infinity in the q plane, it follows that for 
large enough |g|, there is a single leading term |a£(g)| in ( |2.2|) . Hence from eq. (|2.8|) , we see 
that at sufficiently large |g|, the boundary B for the limiting function W{{Kp x G)rb},q) is 
determined by the equality 



IMq-p)] = \Mq-p + 1) 

Note that this is independent of b, so that 



(2.11) 
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B is independent of b for {{Kp x G)rb} 



(2.12) 



(given that the basic condition ( p.5|) is satisfied). Let q = + iqj- One can next enumerate 
the various cases possible for ae{q). The basic theorem that the coefficient of the highest- 
order term, g", in the chromatic polynomial P{Gn, q) of any n-vertex graph Gn is unity 
implies that if a dominant term ai{q) is a polynomial of degree Smax) 

ae{q) = Y.asq' (2.13) 

then 

^ ^max ' ^} 

Consider, for example, the case where ae^{q) is a linear function of q: ai{q) = aiq + ao, which 
reduces to ae{q) = g + ao by ( |2.14| ). Then ( |2.11| ) yields 

«?H=p-(^ + ao) for Smax = I (2.15) 

with q^- undetermined, i.e., a vertical line segment extending to ±ioc + P — (| + «o) in the 
complex q plane. This type of behavior is exemplified by graphs involving Kp adjoined to 
trees or chains of triangles, in which one bond in the Kp subgraph is removed. We shall 
discuss these below. 

If ai{q) is a quadratic function of q, ai{q) = q^ + aiq + ao, then ( p.ll| ) yields an equation 
which has, as its only acceptable solution, 

g^=p-^(l + ai) as \qi\ oo for s^ax = 2 (2.16) 

Hence, the boundary B in this case is, for sufficiently large Ig^l, again a vertical line in the 
complex q plane located at the value of given by ( |2.16| ) and extending to ±ioo. This 
type of behavior is exemplified in our discussion below of graphs involving Kp adjoined to 
chains of squares (i.e., ladder graphs) with various boundary conditions, in which r bonds 
are removed from the Kp subgraph in the manner discussed above. In general, as we shall 
show, if one adjoins Kp to an open chain of /c-gons arranged such that two adjacent fc-gons 
intersect along one of their mutual edges, then the resultant chromatic polynomial has the 
form (|2.8|) with Smax = k — 2 and co(g) = (there should be no confusion in the notation of 



k for the fc-gons and K for Kp). 

This, then, is the algorithm for producing families of graphs depending on three parame- 
ters, p, b, and n, with the property that the limiting function Wr is nonanalytic at 1/g = 0. 
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We have proved this by calculating first the chromatic polynomials for finite graphs and 
then their respective limiting functions W. The key ingredients in the construction are, 
first, the adjoining of the complete graph Kp to Gn, and second, the removal of b of the 
bonds connecting one vertex in Kp to other vertices of Kp. Together, these guarantee, via 
eq. (|2.7|) , that the equation for the degeneracy of the leading term ai is of the form (|2.11| ), 
and the locus of points that solve this degeneracy equation extends to complex infinity in 
the q plane, as we have shown. 

We next give some illustrations of the application of this algorithm. 



III. FAMILIES OF GRAPHS WITH Wr NONANALYTIC AT l/q = 



A. (KpXTn) 



rb 



Perhaps the simplest illustration is for the family Kp x T„ formed by adjoining a succession 
of p vertices to each of the n vertices of a tree graph T„ and to each other. Removing b bonds 
connecting a vertex of Kp to other vertices of Kp in the manner described above yields the 
graph {Kp x Tn)rb- From eq. 



P{{KpXT^U,q) = [f[{q-s)] [{q 



we calculate 

p 
J 

s=0 



P 



\n-l 



b{q-p) 



n-2 



(3.1) 



Our general analysis in eqs. ( |2.11D -( PTTBD applies with 
that, from the general formula ( |2.15| ) we have 



1, ao = —1, and Co(g) = so 



q^=p + - for G = {KpX T)rb 



(3.2) 



Hence, the boundary B consists of the vertical line ( |3.2| ) with —ioo < qj < ioo. The diagram 
describing the regions of analyticity of the limiting function W{{{Kp x T)rb}, q) consists of 
two regions. 



Ri : Re{q) > p + 



(3.3) 



and 



R2 : Re{q) < P + 1^ 



(3.4) 



Mapped to the 1/q plane (a conformal transformation), the image of the vertical line is a 
closed curve which crosses the real axis at the inverse of q^ in ( p.2|) and at the origin. In the 
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1/q plane, the image of region Ri is a compact region enclosed by this closed curve, while 
its complement is the image of the region R2. We find that 



W{{{KpXT)r,},q)=q-v for q E (3.5) 

i.e., Wr{{{Kp X T)rb\,q) = 1 — v/l- For ^ -R2, if q is real, P{{Kp x Tn)rb,q) alternates in 
sign as n increases through even and odd integers, so, strictly speaking, the limit as n — 00 
does not exist. Of course, there is also a corresponding variation in phases in this limit for 
the case of complex q E R2- As we have discussed before, in such a situation, at least the 
magnitude does have a well-defined limit: 

\Wi{iKpXT)rb},q)\ = \q-p-l\ for q E R2 (3.6) 

or equivalently |PVr({(i^p x T)rb},q)\ = |1 — 2±i|. This simple example thus explicitly 
illustrates the nonanalyticity at 1/q = 0; even aside from the choice of the phase in region 
R2, the two expressions above for the magnitude of the reduced Wr functions are different. 
The nonanalyticity in the Wr function thus involves (for any choice of phases in (|1.3| ) in R2) 
a discontinuity in the first derivative dWr/dz, where z = l/qatz = 0. 

Although these are all exact results for the n ^ 00 limit, it is of some interest to see 
how these boundaries develop by studying chromatic zeros for finite n. (Here, we recall the 
definition that the "chromatic zeros" of a graph G are the zeros of the chromatic polynomial 
P{G,q) for this graph.) We carried out these types of studies for a number of families of 
graphs in Refs. |jl3[ and |16|. A general question that we investigated was the following: 



excluding a well-understood subset of chromatic zeros at certain discrete real integer values, 
and considering the remainder, how close are these remaining zeros, for finite n, to the 
boundary B that obtains in the n 00 limit? From our study of the bipyramid graph, in 



Ref. [0] (see also Ref. 0]), it was found that the chromatic zeros near to the real axis in the q 
plane (aside from the discrete zeros at g = and 1 and a zero very near to g = Be^ = 2.618...) 
lie on or near to the arcs forming the boundary curves B{Ri, R3) and B{R2, R3) (defined in 
greater generality in eqs. (|3.11|) , ( p.l2|) below), but the outer zeros do not lie very close to 
the line segments of B{Ri, R2) (given by the p = 2 special case of eq. ( p.l3|) below) and only 
approach these line segments slowly as n increases. We inferred that this latter behavior was 
connected with the fact that this component of the boundary extends to complex infinity 
in the q plane and observed that this type of deviation did not occur for families of graphs 
whose boundaries B were compact and did not extend to complex infinity. 

Here we have carried out an analogous study of the chromatic zeros of {Kp x Tn)rb and we 
find similar behavior as a function of n. Since an interesting feature here is the dependence 
of the locations of chromatic zeros on 6, we focus on this. In Fig. |I| we show the results for 
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p = 4: and n = 18 for b = 1,2, and 3, the full range allowed by (|2.5| ). The outermost zero 
and its complex conjugate do not lie very near to the vertical line with = 9/2 given by 
the p = 4 special case of eq. ( p.2|) . As b increases from 1 to 3, this outermost zero moves 
to larger and slightly smaller g^, and hence slowly toward the above-mentioned vertical 
line. From eq. (|3.1| ), it follows that for p = 4, there are also discrete chromatic zeros at 
g = 0, 1, 2, 3, and 4, and these are evident in Fig. |^. 



8.0 
6.0 
4.0 
2.0 
0.0 
-2.0 
-4.0 
-6.0 
-8.0 



-4.0 0.0 4.0 8.0 12.0 16.0 

Re(q) 



FIG. 1. Chromatic zeros of {Kp x Tn)rb for p = 4, n = 18, and 6 = 1 (•), 2 (+), and 3 (x), where 
the symbols for the points are given in parentheses. In the n — > oo limit, our exact results give 
the continuous region boundary B as the p = 4 special case of eq. (^^), i.e., the vertical line with 
Re{q) = 9/2 and Im{q) arbitrary, independent of b by eq. ( p.l2 ). 
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B. {Kp X Q 



n Irb 



A second illustration is provided by the family of graphs {Kp x Cn)rb obtained by starting 
with the p- wheel Kp x C„ with p > 2, where C„ is the n- vertex circuit graph, and removing 
b bonds connecting a vertex in Kp to other vertices in Kp. From the general formula eq. 
we calculate 



P{{KpXCn)rb,q) = [n(?-^)]x (3-7) 
iq-p+l){q-p~l){{q~p- IT' + (-1)"} + b{q -p){{q~ p^' + (-1)"} 

Our general analysis in eqs. (|2.11|) -( pT5D applies with Smax = 1 and oq = — 1 so that again, 
for sufficiently large B contains the vertical lines extending to ±ioo with i?^ = p + 1/2, 
as was the case with W{{{Kp x T)rb},q)- However, because the Co(g) terms are nonzero in 
this case, the region boundaries in the vicinity of the real axis differ from the simple vertical 



line found for W{{{Kp x T)j,fe},g). By methods similar to those which we used in Ref. flS 
(i.e. working out the conditions for the degeneracy of the leading terms in P), we find that 
the region diagram for W{{{Kp x C)rb},q) consists of three regions: 

Ri : Re{q) > p + and \q - p\ > 1 (3.8) 

R2 : Re{q) < p + ^ and |g-(p + l)|>l (3.9) 

and 

R3:\q-p\<l and |g-(p + l)|<l (3.10) 
The boundaries between these regions are thus the two circular arcs 

B{R,,R,):q = p + e^\-^<e<^ (3.11) 

and 

B{R2, Rs):q = p+l + e'\—<4><— (3.12) 
together with the semi-infinite vertical line segments 

B{Ri,R2) = {q}: Reiq)=p+^ and |Jm(g)| > ^ (3.13) 



10 



These meet at the intersection points 



Qint. 



^ 1 , .73 
p H — ± t — 
^ 2 2 



(3.14) 



The arcs B{Ri, R3) and B{R2, R3) cross the real axis at g = p + 1 and q = p, respectively. 
The region diagram for W{{{Kp x C)rb}, q) is shown in Fig. ||. 



2.0 



1.0 



0.0 



-1.0 



-2.0 



0.0 




FIG. 2. Diagram showing regional boundaries comprising B for W{{{Kp x C)rfe},(?). Breaks in 
the horizontal axis indicate that p is an arbitrary integer > 2. 

Mapped to the 1/q plane, the image of the vertical line is a curve in the right-hand half- 
plane which passes vertically through the origin, enclosing the image of Ri. The enclosed 
region R3 remains enclosed, while the left-hand region R2 in the q plane maps to the region 
exterior to the images of -Ri and -R3 in the 1/q plane. The case p = 2 is the bipyramid family 
of graphs 0, 0. We find that 



W{{{KpX C)rb},q) = q-p for g G i?i 



(3.15) 
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For the other regions, we have, in general. 



\W{{{K.p^CU},q)\ = \q-{p + l)\ for g G i?2 

= 1 for g e /?3 (3.16) 

The Wr functions in regions Ri and R2 are the same as those in the corresponding regions 
for the previous family, {{Kp x T)rb}i while region has no analogue for that family. As 
before, the nonanal5d;icity in Wr is most conveniently discussed in the 1/q plane. There are 
two regions contiguous across the image of B at the origin, 1/q = 0, namely, the images 
under inversion of regions Ri and R2. The nature of the nonanalyticity in Wr at 1/g = is 
the same as that in the {{Kp x T)rb} family. 

Again, it is of interest to calculate the chromatic zeros for some finite graphs in these 
families and see how close they lie to the locus of zeros in the n — 00 limit. We have carried 
out such a study. As an illustration, in Fig. |], we show the chromatic zeros for {Kp x Cn)rb 
with p = 4:, and n = 18 for 6 = 1,2 and 3. This may be compared with the p = 4 special case 
of the plot of the region diagram in Fig. ^ (which is the same for all b, by the result (p.l2 )). 



Note that the discrete chromatic zeros at g = 0, 1,2, and 3 are not part of the continuous 
locus of zeros forming B in the n — » 00 limit. For comparison, the previously mentioned 
bipyramid graph is the case p = 2, for which only the single value 6 = 1 is allowed. We 
find that the complex chromatic zeros near the real axis lie on or close to the arcs B{Ri, R3) 
and B{R2,R3), but the outer zeros do not lie very close to the line segments of B{Ri,R2) 
given by the p = 4 special case of (|3.13|) , namely, =9/2, |g/| > v^/2 and only approach 
these line segments slowly as n increases. Evidently, the behavior of the chromatic zeros as 
functions of b is qualitatively similar to that which we observed in Fig. ||. We have carried 
out analogous studies of the chromatic zeros for other families of graphs constructed via our 
algorithm to have region boundaries B extending to complex infinity and have found similar 
results. 
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4.0 0.0 4.0 8.0 12.0 16.0 

Re(q) 



FIG. 3. Chromatic zeros of {Kp x Cn)rb for p = 4, n = 18, and 6=1 (•), 2 (+), and 3 (x). In 
the n — > oo hmit, our exact results give the continuous region boundary B as the p = 4 special case 
of Fig. |. 

C. {Kp X {Ch)k,n)rb 

We define an open chain of m fc-gons constructed such that a given fc-gon intersects with 
the next fc-gon in the chain along one of their mutual edges as {Ch)k,n where the number of 
vertices is 

n={k-2)m + 2 (3.17) 
The chromatic polynomial for the open chain of m fc-gons is 

PiiCh)k,n,q) = q{q-l)D,{qr (3.18) 
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where 



^.(g) = E(-irf^ (3.19) 



s=0 \ S 



Next, consider the graph Kp x {Ch)^^^ obtained by adjoining Kp to the chain of m fc-gons, 
{Ch)k,n, where n is again given by eq. ( p.l7|) . From (|2.4|) , we have 



P{KpXiCh)k,n,q)= (l[iq-s)]Dkiq-pr (3-20) 

s=0 



The chromatic polynomial P{{Kp x {Ch)k,n)rb,q) is then obtained by substitution of ( p.20|) 
into eq. (p7|) . Note that 

Dk{q) = a,{q) (3.21) 



in eqs. (p78|), ( |2.11|) , and (|2.13| ) general class of graphs. From the above substitution, we 



find that the resultant chromatic polynomial is of the form (|2.8| ) with 

Co(g) = , s^ax = k-2 for {Kp x {Ch)k,n)rb (3.22) 

For sufficiently large \qj\, the region boundary B for W{{{Kp x {Ch)k)rb},q) consists of 
complex-conjugate vertical line segments extending to ±ioo with Re{q) = equal to 

'in=P+ 2{k_2) {{KpX{Ch)k)rb} with A;>3 (3.23) 



where the notation convention of ( |2.1CI| ) is used. Note that if one considers chains of fc-gons 



with k progressively larger and larger, this approaches the limit Yrnvk^ooqn = p + 1/2. As 
defined in the Introduction, the region to the right of this boundary is -Ri, and we denote the 
region to the left as R2- In general (in particular, for > 5) there will be other regions (within 
which W{{{Kp X {Ch)k)rb}, q) is analytic) but we only need -Ri and R2 for our discussion of 
the nonanalyticity at 1/g = since it is only the images under inversion of these two regions 
that border the origin of the 1/q plane. For the resultant limiting functions W we calculate 

W{{{KpX{Ch)k)rb}.q) = Dk{q~p+l)^- for g G i?i (3.24) 
\W{{{KpX{Ch)j,)^rb},q)\ = \Du{q-p)\-^^ for q e R2 (3.25) 

Hence 
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Wr{{{KpX{Ch)k)rb},q) 



1 - 



p — 1 



k-2 

m 

s=0 



-1^ 



{q-p+l)~ 



1 

k-2 



for g G -Ri 
(3.26) 



|lV,({(KpX (C/i)fcM,g)| 



fc-2 



for g G i?2 (3.27) 



For /c > 4, these have the expansions near 1/g = 

k-l 



Wr{{{Kj, X (C/i)fc),J,g) = 1 - p - 1 + 



k-2 



q-^ + 0{q- 



(e.g., 1/g — > 0"*" through real values) and 



\Wri{iK,x{Ch)k)rb},q)\ 



1-ip 



k - 1 

k-2 



q-^ + 0{q- 



as with q E Ri 

(3.28) 

as 1/g — s> with g G -R2 

(3.29) 



These expressions also apply for the case /c = 3 but terminate with the 0{q^^) term, since 
D^{q) is linear and the exponent l/(fc — 2) in eqs. ( |3.24|) and ( 3.25| ) is just unity. Thus, for 
all A; > 3, these results indicate explicitly the nature of the nonanalyticity of the reduced 
function Wr{{{Kp x {Ch)k)rb}^ q) at 1/g = 0. As in our earlier examples, this nonanalyticity 
involves, in general, the sudden onset of a complex phase and apart from this, even for the 
magnitude, a discontinuity in the first derivative dWr/dz at z = 1/g = 0. 

Since the lowest two cases, = 3,4 exhibit particularly simple boundaries i3, we give a 
few explicit results for them. For /c = 3 we find 



P{{K, X {Ch),,nU q) = [n (9 - s)] [{q-p-l){q-p- 2)™ + h{q-p- 1)' 

s=0 



(3.30) 



where from ( 3.17|) , m = n — 2. This has the form of our general eq. (p.8|) with linear a^{q), 
ao = —2 and Co(g) = so that in the n ^ 00 limit, by (|2.15|) , it follows that the boundary 
B consists of the vertical line 



qR=P+2 {(^P ^ {Ch)3)rb} 

The general formulas ( p.26| ) and ( p.27| ) reduce to the simple expressions 



Wr{{{Kj,x{Chh)rb},q) = l 



P+1 



for q E Ri 



(3.31) 



(3.32) 
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and 



\Wr{{{Kp^{Ch);)r,},q)\ 



p + 2 



for g G -R2 



(3.33) 



For the = 4 case, D/^i^q) = ai{q) = — 3g + 3, so that ( p.23|) or the quadratic special 
case eq. (|2.16|) with ai = —3 and co(g) = apphes and yields 



q,=p + l for Wi{iKpXiCh),U},q) 



(3.34) 



Using ( |3.20| ) and ( p.7| ), we determine that the region boundary of W{{{Kp x ((7/1)4)^6}, g) 
is again precisely the vertical line with given by ( |3.34| ) and Ig^l arbitrary. Note that this 
holds even though Di{q) = a^{q) is quadratic. However, the generic behavior for higher k 
is that in the vicinity of the real axis, B is more complicated, and the complex-conjugate 
vertical line segments with q^ given by (|3.23|) apply for Ig^l > k^, where is a fc-dependent 
constant. 

From eqs. ( |3.26| ) and (|3.27|) we have 



Wr{{{K, X {Ch)^U}, q) = [l- (2p + l)q-' + {p'+p+ l)q 
and 



1/2 



for qeRi (3.35) 



1 - (2p + 3)g-i + (p2 + 3j9 + 3)g- 



1/2 



for g G i?2 (3.36) 



D. {Kp X Ln^bcjrb 

A slightly more complicated illustration is provided by starting with an n-vertex ladder 
graph, i.e., chain of squares, as in the k = A case discussed above, but with periodic or 
twisted boundary conditions (rather than open boundary conditions) denoted Ln,pbc and 
Ln,tbc, respectively. Note that any even number of twists is equivalent to no twist and any 
odd number of twists is equivalent to one twist. Following our general algorithm in ( |2.7| ), we 
adjoin Kp to this ladder graph and then remove b of the bonds connecting a vertex in the 
Kp graph to other vertices in Kp. Using the basic results 

P{Ln,pbc, g) = (g' - 3g + Sr/' + (g - 1){(3 - g)"/2 + (1 _ g)"/2| + ^2 _ 3^ ^ ^ ^3 37^ 

and 

P(L„,,„ g) = (g2 - 3g + 3)"/^ + (g - 1){(3 - g)"/^ - (1 - g)"/^} - 1 (3.38) 



16 



we can apply our general analysis above. First, we note that because the terms raised to the 
n'th power, which can be leading terms as discussed in conjunction with eq. (p.2|) , are the 
same for P{Ln,pbc,<l) and P{Ln,tbc,(l), it follows that 

WiiiKp X LpbcU}, q) = Wi{{Kp x LtbcU}, q) (3.39) 

and hence the boundaries B are identical for these two functions. In both cases, from our 
general results, it follows that B contains complex-conjugate vertical line segments that 
extend to ±icx3 with 



q^=p+l (3.40) 

as for the case of free boundary conditions, eq. (p.34|) . These line segments extend outward 
to qi = ±00 from the intersection points 

qint., q*int. = P + 1 ± iVl (3.41) 

At these intersection points, the boundary bifurcates into two, which extend down and cross 
the real axis at 



qR,cross Vl P 



+ V2 (3.42) 



In Fig. ^ we show a plot of the region diagram for W{{{Kp x Lpbc)rb},q) = W{{{Kp x 
Ltbc)rb},q)] again, by ( |2.12|) , it is independent of b. 
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FIG. 4. Diagram showing regional boundaries comprising B for W{{{Kp x Lxbcjrb}^ q) where xbc 
denotes periodic or twisted boundary conditions {pbc, the). Breaks in the horizontal axis indicate 
that p is an arbitrary integer > 2. 

IV. A GENERAL CONDITION GOVERNING THE (NON)ANALYTICITY OF 

Wr{{G},q) AT l/g = 

As we stated at the end of section 0, the key ingredient in our algorithm to construct 
families of graphs {G} with the property that B extends to complex infinity in the q plane 
and hence that WT{{G},q) is nonanalytic at 1/g = is to produce a chromatic polynomial 
with the feature that there are two leading terms with a degeneracy condition ( |2.11| ). From 
(p.7|), these leading terms are actually the same function, just evaluated at arguments that 
differ by unity, as is evident in eq. ( p. 111) . Here we generalize from this study to state the 
following theorem: 
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Theorem. Consider a family of n-vertex graphs and its n — oo hmit {H}. With P{Hn, q) 
in the form of eq. (|2.2| ), the region boundary B of the resultant asymptotic limiting 
function W{{H},q) extends to complex infinity in the q plane, and hence Wr{{H}^q) is 
nonanalytic at 1/g = in the 1/q plane, if and only if the locus of solutions of the 
degeneracy condition of leading terms 

\aM)\ = Wi{q)\ (4.1) 

extends to complex infinity. 

Of course, this locus of points obeys the condition ( |2.1|) . Our algorithm yields the family of 
graphs {Kp x G„)rb, depending on p and 6, whose n ^ oo limit, {{Kp x G)rb} satisfies the 
above condition. 



A general feature of our results is that for the families of graphs we have constructed and 
studied, for which B extends to complex infinity in the q plane, the image, under inversion, 
of B passes through the origin of the 1/q plane with an infinite tangent, i.e., vertically. This 
reflects the property that the portion of B that extends to complex infinity is comprised of 
a vertical line segment and its complex conjugate, with a fixed value of Re{q). 

One salient feature of our study is clearly that none of our families of graphs with 
Wr{{G},q) nonanalytic at 1/q = is a regular lattice graph. Our results are therefore 
consistent with the assumption underlying the original series calculations, that a sufficient 
condition for Wr{{G}, q) to be analytic at 1/g = is that {G} be a regular lattice graph A. 
We state this formally as the following conjecture: 



Conjecture: Let {G} denote the infinite-n limit of a family of graphs Gn- A sufficient 
condition for the resultant reduced functions Wr{{G}, q) to be analytic at 1/g = is that 
{G} is a regular lattice graph {G} = A. 



It is clear from Ref. that the property that {G} be a regular lattice is not a necessary 
condition for the associated Wr{{G}, q) to be analytic at 1/g = 0; in that work we calculated 
iy({G}, g) functions for a number of families {G} which are not regular lattice graphs, but 
for which the corresponding Wr{{G},q) functions are analytic at 1/g = 0. 

For a regular lattice A with coordination number a natural reduced function is defined 

by 

iy({G},g)=g(l-g-ipW(A,i/) (4.2) 
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for which the large-g Taylor series expansion is |T0|- ]T2 



CO 1 

W(A, y) = 1 + E «^A,n2/" , y = —— (4.3) 

n=l 1 ^ 

Clearly, W{K,y) is analytic at ?/ = if and only if Wr{Is^,q) is analytic at 1/g = 0. In 
Refs. Ijl^, [0, [|18| we have made detailed comparisons of existing large-g series expansions 
with high-precision Monte Carlo measurements of W^(A, q) as well as rigorous lower bounds 



that we have derived and have found excellent agreement for g > 4 on a number of 
different lattices, including both homopolygonal (e.g. square and honeycomb) lattices and 
heteropolygonal Archimedean lattices (composed of regular polygons of more than one type, 
such that all vertices are equivalent), viz., the 4 ■ 8^ lattice, for which we have calculated 
a large-g series). This excellent agreement provides motivation to include lattices involving 
packings with different regular polytopes under the term "regular lattice" in the above 
conjecture (here, polytope is the general mathematical object that subsumes the polygon in 
2D and polyhedron in 3D pi). 



V. CONCLUSIONS 

In this paper we have addressed a fundamental problem in graph theory with important 
implications for statistical mechanics, namely the question of the analyticity of Wr{{G},q) 
at 1/g = 0. In order to understand the phenomenon of nonanalyticity of this function 
at 1/g = better, we have constructed a general algorithm for producing infinitely many 
families of graphs, each depending on two parameters p and 6, with Wr functions that are 
nonanalytic at 1/q = 0. We have studied the properties of several of these families. We 
have also stated a general necessary and sufficient condition on the chromatic polynomial 
of a family of graphs such that the resultant Wr{{G},q) is nonanalytic at 1/q = 0. This 
condition explains the source of the nonanalyticity in the cases where it occurs. The results 
of our study are consistent with the conjecture that a sufficient condition (we know that this 
is not a necessary condition) for Wr{{G}, q) to be analytic at 1/g = is that {G} = A is a 
regular lattice graph. 

This research was supported in part by the NSF grant PHY-93-09888. 
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